home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / IMTQL1.z / IMTQL1
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. ____IIIIMMMMTTTTQQQQLLLL1111((((3333FFFF))))                                                        ____IIIIMMMMTTTTQQQQLLLL1111((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      IMTQL1, SIMTQL1  -  EISPACK routine.  This subroutine finds the
  10.      eigenvalues of a SYMMETRIC TRIDIAGONAL matrix by the implicit QL method.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee  iiiimmmmttttqqqqllll1111((((nnnn,,,, dddd,,,, eeee,,,, iiiieeeerrrrrrrr))))
  15.           iiiinnnntttteeeeggggeeeerrrr          nnnn,,,, iiiieeeerrrrrrrr
  16.           ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn dddd((((nnnn)))),,,, eeee((((nnnn))))
  17.  
  18.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssiiiimmmmttttqqqqllll1111((((nnnn,,,, dddd,,,, eeee,,,, iiiieeeerrrrrrrr))))
  19.           iiiinnnntttteeeeggggeeeerrrr          nnnn,,,, iiiieeeerrrrrrrr
  20.           rrrreeeeaaaallll             dddd((((nnnn)))),,,, eeee((((nnnn))))
  21.  
  22.  
  23.  
  24. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  25.      On INPUT
  26.  
  27.      NNNN is the order of the matrix.
  28.  
  29.      DDDD contains the diagonal elements of the input matrix.
  30.  
  31.      EEEE contains the subdiagonal elements of the input matrix in its last N-1
  32.      positions.  E(1) is arbitrary.  On OUTPUT
  33.  
  34.      DDDD contains the eigenvalues in ascending order.  If an error exit is made,
  35.      the eigenvalues are correct and ordered for indices 1,2,...IERR-1, but
  36.      may not be the smallest eigenvalues.
  37.  
  38.      EEEE has been destroyed.
  39.  
  40.      IIIIEEEERRRRRRRR is set to ZERO       for normal return, J          if the J-th
  41.      eigenvalue has not been
  42.         determined after 30 iterations.  Calls PYTHAG(A,B) for sqrt(A**2 +
  43.      B**2).  Questions and comments should be directed to B. S. Garbow,
  44.      APPLIED MATHEMATICS DIVISION, ARGONNE NATIONAL LABORATORY
  45.  
  46.  
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.